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Abstract: Disformal transformations have proven to be very useful to devise models of 
the dark sector. In the present paper we apply such transformation to a single scalar 
field theory as a way to drive the field into a slow roll phase. The canonical scalar field 
Lagrangian, when coupled to a disformal metric, turns out to have relations to bimetric 
dark matter theories and to describe many specific dark energy models at various limits, 
thus providing a surprisingly simple parametrisation of a wide variety of models including 
tachyon, Chaplygin gas, K-essence and dilatonic ghost condensate. We investigate the 
evolution of the background and linear perturbations in disformal quintessence in order to 
perform a full comparison of the predictions with the cosmological data. The dynamics 
of the expansion, in particular the mechanism of the transition to accelerating phase, is 
described in detail. We then study the effects of disformal quintessence on cosmic microwave 
background (CMB) anisotropics and large scale structures (LSS). A likelihood analysis 
using the latest data on wide-ranging SNIa, CMB and LSS observations is performed 
allowing variations in six cosmological parameters and the two parameters specifying the 
model. We find that while a large region of parameter space remains compatible with 
observations, models featuring either too much early dark energy or too slow transition to 
acceleration are ruled out. 
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1. Introduction 

Observations of the Cosmic Microwave Background temperature anisotropy reveal that a 
mysterious constituent with negative pressure, so called dark energy, accounts for about 
70% percent of today's mass energy budget and is causing the expansion of the universe to 
accelerate [|l|, These observations are in remarkable concordance with the observations 
of distant supernovae ^. A major challenge in present day cosmology is to discover 
the physical nature of dark energy. Meanwhile, though the evidence for the existence of 
dark matter has been accumulating for several decades, the quest to find its precise nature, 
whether as a particle described by an extension of the standard model, a more exotic field or 
even a modification of General Relativity, is ongoing. Many possibilities have been explored 
in attempts to explain the workings of the dark sector ||5|, ^, ^, ^, ^, |l^, ll|, 1^, 13, 14, 15|, p^. 



The phenomenologically simplest but theoretically very problematical dark energy is 
the Einsteins A term [^]. Quintessence provides a dynamical alternative to the static 



cosmological constant |T8| , 19 j. It tracks or scales with the background energy density. 



therefore more naturally resulting in similar orders of magnitude for the dark energy and 
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dust-like matter energy densities today. However, one needs some mechanism to end this 
scahng and to onset the acceleration. Several possibilities have been considered, in par- 
ticular introducing a suitable bump into the form potential [^], taking into account the 
Gauss-Bonnet term [21, B^, coupling the field to other matter 24 1, considering non- 



canonical Lagrangian |25|, or introducing many fields [26|. In the following, a disformal 



relation |27] is applied for this purpose. 

Disformal relations have recurrently appeared in cosmology in models of alternative 
dark matter, in particular TeVes p^] and other relativistic and covariant formulations of 



MOND [^], and in bigravity theories [30|. This comes about for several reasons: theoreti- 
cally, it is related to the Born-Infeld type of Lagrangians (as will become clear below), and 
phenomenologically, it is needed to produce the observed lensing without resorting to par- 
ticle dark matter [^] . Let us mention that the possible relevance to inflation has also been 
considered [^], and that disformality is a key to varying speed of light, which has been 
considered as an alternative to inflation [pSl p41 . The usefulness of conformal transforma- 



tions is well known and appreciated |35, 36, 37 1, but many aspects of the wider framework 



remain uncovered. The conformal transformation relates the simplest scalar-tensor theo- 
ries (including the f{R) theories |^8|) to general relativity ||39[| , but in the general case one 
is prompted to turn into extended classes of transformations. 

Here we make a preliminary step along this direction by exploring the perhaps simplest 
possible set-up: the canonical scalar field coupled to a disformal metric. It turns that 
already then many relations between seemingly disconnected models emerge. If the scalar 
field is reduced to a constant term, it appears as a tachyon or a Chaplygin gas in 
the physical frame. When the kinetic term is included, some other previously considered 
models can be recovered, as we will discuss in Section |2|. In the present study we then 
focus on one simple possibility, disformal quintessence. We point out that it can act as 
viable dark energy when the parameters in the Lagrangian are of the Planck scale. Section 
m is devoted to study of this model, and contains the main results of this paper. We 



describe in detail the background evolution of this model in subsection 3T. Particular 
attention is given to the transition mechanism providing an exit from the scaling era. The 
details of this transition depend on the two parameters of the model, and thus they can 
be constrained by the SNela data. We also consider the evolution of perturbations in 



subsection 3.2 in order to compute the CMB and matter power spectra. Armed with these 



solutions, we perform a Monte Carlo analysis of the model in subsection 3.2, combining 
data from many different cosmological experiments. Wide parameter range is found to 
be compatible with observations, but certain parameter region, corresponding to shallow 
slope of the exponential potential or of the disformal factor, can be ruled out. We conclude 
briefly in section ^ 

2. Dark energy and unified models from disformal relation 

A generic scalar field Lagrangian is a function of two variables, 

£ = £(<A,X). (2.1) 
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Expanded as a series in the kinetic term X = ^g^^ it reads 

C = C{4>, 0) + 0)X + C,xx{(k, 0)^2 + . . . (2.2) 

The first term corresponds to the potential, and the second term, if canonically normal- 
ized, equals one. The dimension of the prefactor in third term is and the relative 
importance of this term seems to be H/M. Unless the the prefactor is huge, this term is 
negligible in the late universe. Apparently a canonical scalar field should be an adequate 
description of almost any natural case. However, we will soon see that this can be avoided 
for suitable forms of C. This comes about simply because we can have an exponential 
amplification of the coefficients in the formal series ( |2.2| ). 

We consider the following generalization of the conformal transformation 

Qap = A{4>)9o,p + B{cl))4>^a4',p (2.3) 



defined by a scalar field (f>. Equation (|2.3| ) constitutes the most general relation between 
metrics that preserves causality and covariance and is defined through (p and its first deriva- 
tives only [27|. 



2.1 Disformal dark energy 

The simplest theory where ( |2.3| ) appears to be relevant is 



Sa = - J dJ'x^/^K. (2.4) 

It turns out that such an action corresponds to a matter component obeying the equation 
of state 

P = • (2.5) 

P 

Thus the cosmological term of a disformal metric is a tachyon in a constant potential^ . This 
describes unstable branes in string theory The equation of state (2.5) corresponds 



to Chaplygin gas , which has been applied in cosmology in attempts to unify dark matter 
and dark energy [^]. Our result shows that this rather exotic fiuid can be obtained in an 
extremely simple way from a cosmological constant term by taking into account its possible 
disformal relation to the matter metric. 

The A term describes the potential of the tachyon, and therefore it is natural to consider 
it as a field dependent term. A straightforward generalization of this is to include also a 
kinetic term for the field. In the case of canonical kinetic term we obtain the model 



-g^^g^'cp^^^^u + Vi^)). (2.6) 



As it does not complicate our analysis we keep the sign of the kinetic term, e arbitrary. 
The case e = 1 was considered in Esl as a variation of the quintessence scenario. The effect 



^The action given by (2.4) can be easily written in terms of the gravitational metric using (2.17) for the 



determinant of the disformal metric. 
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of the disformal relation is then to end the scahng era and begin the accelerated expansion. 
The canonical field in the disformal metric assumes the rather unappealing form in terms 
of the physical metric 



For derivation of this form and a variation of the model, see subsection 2.4. Then, expansion 
in the kinetic energy follows 

£^ = + AX{£ + BV) - ^BX^{2e + BV) (2.8) 

+ ^(3e + W) + ^(2e + W) + ... 
Due to the form of the expansion all terms after first order in X are of the form B'^~^ A'^~'^ 

IpnE + QnBV] with 

PmQn rational numbers. The argument following ( p.2| ) can be circum- 
vented if the coefficients are adequately enhanced so that the expansion is not valid any- 
more. One possibility is to choose an exponential form for the disformal function B (with 
A = l). 



Notice that the relation (2^) becomes meaningful only after we have specified the 
coupling to gravity. We consider the gravity sector to be given by the Einstein-Hilbert 
term of the unbarred metric. Then action ( |2.6D could be considered to be the Einstein 



frame version of a scalar tensor theory with derivative couplings |44]. In terms of the 
Jordan frame metric, the variation of the action ( |2.(]| ) alone gives the standard equation of 
motion of the field, □(/) + V'{(f)) = 0, and the canonical stress tensor - however, one would 
have to take into account that the gravity piece would then be nonminimally coupled to 
the scalar field, resulting considerably more involved equations. 

As we have seen, the disformal relation has a very simple connection with the Chaplygin 
gas. It turns out that the theory given by ( p.6| , p^ ) can encompass a wide range of dark 
energy models through different choices of e, B, V and A. Some variations are shown in 
Table |l[ Note that the very inclusion of the disformal metric g^j^i, makes the scalar field 
dynamical, as it contains derivatives of 0. Disformal quintessence will be discussed in depth 
in section |^. 

2.2 Bigravity 

Recently the Born-Infeld form of the action has been considered in the framework of bi- 



gravity theories |30|. Following Ref.|51|, let us write 

S = Seh{9) + CSEHiq) + j V^q<Px {c^ - ^Aql^'g^u^ , (2.9) 

where the two first terms are the Einstein-Hilbert actions of the metric g and the metric 
q, respectively, with a relative strength given by C. The interaction term coupled to the 
metric q includes an additional scalar field as 

£^ = ^Bq>'''4>^^4>^, + U{ct>). (2.10) 
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S((/>) Model 



Remarks 



A Cosmological constant w = 

1 Vi^P) Quintessence w > 
-1 y{4>) Phantom quintessence w < 



-1 

-1, lll,|9| 
-1 






> 


±A 


(Anti-)Chaplygin gas 





1 


v{4>) 


Tachyon condensate 


-1 


^(0) 





K-essence ^ 


-1 







Dilatonic ghost ^ 


-1 






Disformal phantom 



Eq. H @ 



mm 

^ = S > 0, H, |5| 

mm 



1 > y{4>) Disformal quintessence |43|, Section 



Table 1: Disformal dark energy models obtained from the action ( |2.6| ). The first three models do 
not use disformal transformation. The second set contains examples which have been thoroughly 
studied in the literature. Disformal quintessence will be considered further in this work, while the 
phantom version is left for further studies. All models are purely disformal {A = 1) unless explicitly 



mentioned. A modification to obtain further variations is briefly discussed in subsection 2.4 



In general, A, B and C could be functions of though we do not write this explicitly. By 
varying with respect to q^^'^ , one finds the relation of the determinants 



(2.11) 



where K^j^^ is the Ricci tensor and K the Ricci scalar associated to g^j^ and g is the deter- 
minant of the metric 



Using the trace of the field equations for q^^ we further get that 

2^ ' 



n 



2nU 



9- 



(2.12) 



(2.13) 



This leads us to the observation that the scalar-bitensor theory (^) we started with is 
equivalent to 



S = Seh{9) + / d^x 



n 



2nU 



(2.14) 



Setting A = U = 1, B = we get the Einstein-Born- Infeld (EBI) theory. This slightly 



generalizes the derivation presented in |51], where the equivalence of a bimetric theory and 
the EBI theory was pointed out. It is interesting that the disformal scalar field enters 
here into play, since up to now the EBI has been shown to mimic dark matter both at 
galactic and cosmological scales, but the structure formation having tension with the data 
when simultaneously providing dark energy The Chaplygin gas ensuing from the 

tachyon Lagrangian on the other hand can act as dark energy but the structure formation 
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bounds rule it out as a candidate for dark matter |53, 54 1. This suggests that the successful 
unification of the dark sector ultimately might require considering a theory along the lines 
of ( |2.9| ) . In the present study, we make a step towards this by analyzing in detail the scalar 



sector of such a theory, which could be responsible for dark energy. 
2.3 Rapid transition 

We also note that a solution to the structure formation problem of the unified models 
discussed above has been suggested based on a rapid transition of the unified fluid from a 
dark matter-like phase to an effective cosmological constant |^. Then, the sound speed 
could be nearly vanishing except for the transition period that is very short in cosmological 



time scales. As we will show in detail in section 3.1.3, the disformal relation can naturally 
provide means for such a sharp transition, thus suggesting that a simple field theory may 
generate the desired scenarios. However, to assess their viability, one would have to consider 
the full nonlinear evolution of cosmological structures, since the linear perturbation theory 
predicts that when the cosmological constant phase is reached, the structures that have 
formed earlier are simply wiped away, resulting in vanishing matter power today. In reality 
though the already virialized objects may be expected to retain their form as the field 
changes its nature at the very largest cosmological scales, but the details of this are beyond 
the scope of our present study. 

2.4 A variation of the model 

By contraction, the inverse disformal metric can be shown to be given by 

r' = \( - -ir^i^r^^'') ■ (2-15) 

Then we see that 



AV A + 2BX 



X = —^ ■ (2.16) 

A + 2BX ^ ' 

The determinant of the disformal metric is (for a method to obtain the result see Appendix 
C of Ref.||]) 



--g = A^sj\-r2^X^g. (2.17) 



Using these formulas the form (2.7) follows immediately. 

One may also consider an alternative prescription for the canonical field (^]^). The 
kinetic term is written there is terms of the inverse metric as X = ^g'^'^cft^ij^cl)^,^. The 
alternative formulation employs the metric in the combination X = ^g^u4''^4'''^ ■ This is a 
consistent but not the minimal prescription, since depending on the viewpoint, we are doing 



one of the following: 1) mixing the two metrics (since the derivative indices in ( ^.17 ) are 



raised with the unbarred metric) 2) considering non-canonical field (the unbarred metric 
can be barred by introducing field combinations) 3) defining the kinetic energy in terms 
of differentials with respect to one- forms (and not the coordinate vector as usually). Now 
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X = X{1 + 2BX), and in terms of the matter metric, this Lagrangian corresponding to 
this model is 

= A^VA + 2BX [X{1 + 2BX) + V{(t))] . (2.18) 
We shall not consider this variation further here. 



3. Disformal quintessence 

We will now explore the model given by action (^]^) and ( |2.7D with e = 1. For the sake of 
simplicity, we will restrict to the purely disformal case A = 1. The dynamical equations 
derived in and hold for any choice of B and V, but we will rather focus on an 
exponential dependence for both functions 

V = M^ e-^'i'l^v . (3.2) 

The model parameters are chosen such that the disformal features are negligible at early 
times so the field behaves as normal quintessence. The form for V ensures the existence 



of scaling solutions in this regime if a > 2 [57]. The exponential form for B allows the 



VC\ 2£ 



disformal features to become relevant without introducing a new scale. 
3.1 Background evolution 

For a flat FRW metric, the energy and density pressure in the Einstein frame read 

,= '(.^ + v) (3.3) 



where the lapse function 



measures the relative time flow in the disformal frame ( |2.3| ) relative to the Einstein frame 
i-e. 900 /goo- The disformal factor 

P = B(/>2 , (3.6) 

is a measurement of the deviation with respect to canonical quintessence. 

The nonlinear field equation can be written in a form analogous to a single harmonic 
oscillator with coefficients that depend on (j), (j) 



P = VC(e^-V) (3.4) 



£ = 1 - 5(^2 , (3.5) 
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The analogs of mass, friction and potential terms can be written in terms of the disformal 
factor and the lapse as 

M = e{l + ^V)+£BV , (3.8) 

T = C[eil - ^V) + CBV] , (3.9) 

V = [e^<P^ + ^CV]B,^<p^ + . (3.10) 

In addition to the disformal factor and the lapse, we find the dimensionless disformal- 
potential factor BV as yet another measurement of the disformal properties. 

3.1.1 Tracking 

We assume that the field has reached the tracking solution with VL^ = njc? before any 
relevant time. This behavior will last while T) and BV ^ 1. If the field is tracking the 
matter density and I? 0, then Y ^ }^(\?' and both conditions hold simultaneously. The 
disformal factor evolves as 

= 2^ [S^gpMp-^] oc 0^(1-^) , (3.12) 

and will be a growing function of time only if /3/a > 1, a condition necessary to push the 
field out of the attractor. The parameter produces a shift in the transition time and 
hence fi,^ is a monotonous growing function of it. In the following we will discuss the model 
in terms of a, /3/a and 

3.1.2 Accelerated expansion 

When the disformal factor reaches 0(1) values, the lapse function reduces significantly and 
the field, as seen in the Einstein frame, is pushed towards a slow roll phase. This produces 
a negative equation of state responsible for the acceleration of the universe. Figure ^ 
displays the effects of the transition in for different values of /3/a. High values lead 
to a rapid slow down of the field and a more negative equation of state while low values 
are associated with gradual and longer transitions in which the field preserves a significant 
velocity at later times. The reasons for this behavior and the equation of state will be 
analyzed in further detail below ( |3.1.3 ). 

The disformal factor depends on both 0, and has a nontrivial behavior (See Figure 
|2|). The dynamics are similar in all cases, but become more clear for higher values of /3/a. 
The transition is associated with a maximum and a rapid fall towards a certain value, 
followed by a slower reduction. The existence of a maximum follows from the role of T) in 
the disformal metric, which implies that 2? < 1 at any time. Otherwise the metric would 
run into a singularity and many quantities would blow up, notoriously the energy density 
( |3.3| ). A slow down of the field is therefore to be expected whenever T> approaches unity. 

No similar bound can be put on the potential disformal factor, which grows exponen- 
tially with (/) and can overcome the small velocity and reach large values after the freeze 
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Figure 1: Disformal quintessence dependence /3/a as a function of redshift [a — 10 and fiq — 0.7). 
Left: Equation of state for the field. Higher values of /3/a produce sharper transitions while low 
values lead to smoother ones.. Right: Speed of sound squared (see section 5.2), equivalent to 
TjM « 1 - 2?. 




Figure 2: Disformal functions P, BV for different /3/a as a function of the scale factor (a = 10 and 
= 0.7). Before the transition both lines evolve together. The disformal factor V < 1 corresponds 
to the lower lines. The bounce avoids a metric singularity in (2.3). The disformal-potential term 
BV corresponds to the lines that cross above one. This term dominates 7W after the transition. 
Note that the evolution has been extrapolated to a = 10. 



out. It can be seen how the change of slope in BV roughly corresponds to the transition 
(i.e. the slow down of the field). This factor is the dominant contribution to M and 
will produce a large inertia for the field. Despite the disformal factor is restored after the 
transition, the dynamics of the field might be strongly suppressed by the large value of 
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BV. 



3.1.3 Disformal transition in detail 

If f3 <^ a, the disformal factor is well below one except during some interval around freeze 
out. This observation suggest an expansion to first order in D to gain some intuition about 
the model. Let us regard the terms in (|3.7| ) in the light of this expansion. For the friction 
term we get a simple correction with respect to the usual quintessence case [T j M. = 1) 

^ = \-V^O{p'). (3.13) 

The expansion olV j ^A taking into account the exponential dependence on the field 
yields 



V a V 



M Mpl + BV 



1 - 2BV 

1 + 3P- 



2(1 + BV) 



R 1 iX'^ _L y 



The first term reduces to the usual V^^ at early times, but it will receive an additional 
suppression by the BV factor in the denominator growing large at later times (The ex- 
pression in brackets remains of order 1). The second term represents the purely disformal 
contribution to the potential and is equally suppressed by BV at late times. Note that it 
has the proper sign to cause a deceleration of the scalar and is suppressed by the disformal 
factor prior to the transition. 

Neglecting fine details from ( p. 14 ), the force acting on the field will flip sign when 



(3.15) 

This force associated with the transition can grow much larger than V^^ if /3 ^ a, which 
is the ultimate cause for the transition to begin so rapidly. As the field slows down, the 
disformal factor oc (fP' and the restoring force are reduced. Eventually, the usual term will 
come to dominate the force again when ( p^ is fulfilled again. We can see in the plots 
that the disformal factor stabilizes approximately around this value in Figure 

As the transition happens the field increases its value and the BV factor in the de- 
nominator strongly suppresses the force term acting on the field. The later and softer 
slow down would then be driven by the friction term, which does not suffer this or other 
suppression. 

3.1.4 Equation of state 

We can compute the equation of state from equations (0, 3^) to first order in (jP' /V and 
T> to he w —\ + T> + ^(jP /V . At late times, when the kinetic energy is negligible, the 
equation of state becomes 

w^-l + V, (3.16) 

and we expect that wq + 1 ^ a/ (5 according to the discussion above. Figure |3| confirms this 
trend for a wide range of parameters. For larger amounts of early dark energy the freeze 
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1+Wn 



0.1 



0.01 



0.001 




Figure 3: Field equation of state at present time as a function of (3 /a for different values of a. 
There is a clear trend towards — f + a/ (3 unless a or /3/a are very low. For the red curve the 
amount of early dark energy is very large and the disformal transition is not finished at present 
time, hence the deviation. 



out does not finish before the required is reached, as it becomes clear for the curve with 
a = 3. Similar considerations apply for long lasting freeze out (low values of /3/a). 

The properties of the accelerated expansion are therefore linked to the quotient /3/a. 
Apparently the only role of a is to regulate the amount of early dark energy and therefore 
shift the transition time to compensate the time interval necessary to achieve ^a, today. 



3.2 Perturbations 

We considered the perturbations in the synchronous gauge [^8| . The relevant equations are 
the same as in were p is given by the Lagrange density ( ^.71) . With the definitions from 
the previous section, we can write the energy density and velocity perturbation induced by 
the disformal field as 



^-5/2 



M(l)5(l) + V5(t) 



The Klein-Gordon equation obeyed by such perturbations can be written as 



+ 



3i/ + 



M 
M 



+ 



M 



1 ; : 



0. 



(3.17) 
(3.18) 

(3.19) 



hk is the interaction with gravity introduced through the covariant derivative. The field 
mass generalizes the term and is given in terms of (background valuated) partial 



derivatives of (2.7) as 



1 



(3.20) 
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3.2.1 Tracking 



As in the homogeneous hmit, the smallness of the disformal factor ensures no significant 
departure from the canonical quintessence case at early times. 

The speed of sound characterizes the propagation of fluctuations within the field. For 



the purely disformal case it is given by [43| 



The allowed values of BV and T) bound its value between and 1 (See Figure ||). It is very 
close to one before the transition, according to the canonical quintessence limit. Therefore, 
the growth of perturbations at the tracking stage will be strongly suppressed and by the 
time of the transition the dark energy overdensity will be many orders of magnitude below 
the matter contrast. 

As for usual quintessence, the finite amount of early dark energy will increase the 
expansion rate and damp the formation of structure in all epochs. This departure of 
matter domination is refiected in a slope in the Newtonian gravitational potentials and a 



slight enhancement of the first acoustic peaks on the CMB. pH, 61]. 



kiQ = 1 00 kiQ = 1 




Figure 4: Large scale perturbations for a = 10. Top: Dark energy power spectrum (arbitrary 
units). Lower values of {3/ a allow more clustering due to different expansion rate and speed of 
sound. Bottom: Gauge invariant potentials for the same fc— modes. The main differences with A 
are departure from matter domination due to early dark energy and the different times in which 
the acceleration epoch begins. 
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3.2.2 Accelerated expansion 



The transition causes changes in the terms involved in ( 3.17 , 3-18 , 3.19| ) such as the sign flip 
and rapid growth in V associated with the freeze out. This variations are compensated 
by corresponding changes in 5(j), which reduces in magnitude and velocity to maintain 
the energy density at a similar level (see Figure Q). The speed of sound drops below 1 
for a certain time and then stabilizes with a value that depends on /3/a (See Figure |l]). 
This would in principle enhance the clustering properties of the disformal field after the 
transition, but the inhomogeneities in the fluid can grow no faster than the matter ones, 
and since they are many orders of magnitude smaller at the beginning of the transition 
they do not play a relevant role. Besides, the transition is associated with the accelerated 
expansion and therefore the growth of structure is further suppressed. 

The main departures from ACDM are caused by the different expansion rates in the 
accelerated phase and the different transition times. Lower values of render the transi- 
tion slower, but a longer expansion period is required to achieve the same density fraction. 
This both affects the Integrated Sachs Wolfe effect and shifts the angular scales of the CMB. 
Higher values mimic the A behavior at late times, and will only reflect the departures from 
matter domination described above for small a, i.e. significant early dark energy. 



a= 10 




Figure 5: CMB spectra for different values of (3 /a (a = 10) and relative deviation with respect to 
A. The initial amplitude of the perturbations has been chosen to correct the amount of early dark 
energy and fit the first peak. The departures from ACDM are noticeable as a shift of the angular 
scale and the increase of the Sachs- Wolfe plateau, but is only significant for very low values of /3/a. 
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a= 10 




0.0001 0.001 0.01 0.1 1 
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Figure 6: CDM spectra for different values oi (3/a {a = 10) and relative deviation with respect to 
A. The spectra fail to converge to A in the limit of high (3 /a due to the finite amount of early dark 
energy at a = 10, which damps cosmic growth. As the shape of the spectrum remains the same, it 
is difhcult to obtain constraints from it due to the ignorance of the linear galaxy bias. 



3.3 Constraints 

A Monte Carlo Markov Chain simulation was used to obtain constraints on the parameter 
space of the theory by comparing with WMAP 7 year data [0, supernovae from the Union 
dataset SDSS DR7 baryon acoustic peak position [^] and matter power from the SDSS 
luminous red galaxies sample [^. The uncertainty in the linear galaxy bias factor was 
treated by choosing the highest likelihood value in the range (1,3). This choice accounts for 
the uncertainty associated while avoiding too large bias factors. Supernova data included 
systematic errors and possible spatial curvature was neglected. The upper bounds in 
the model parameter space a G (3,20) and (3 /a E (1.5,30) was chosen by sensitivity 
considerations, since very little departure from A can be observed at these values. The lower 
bounds ensure the existence of matter attractor solution for the field {a > 3), disformal 
transition (/3 > a) and avoid computational difficulties. The simulation was performed 
through a modification of the publicly available Boltzmann code CMBeasy |6^. Four 
chains were run adding up to 20870 accepted models. An additional MCMC using only 
background information (SNe, BAO and CMB distance priors as described in section 5 of 
Q) was run for comparison. Results are displayed in Table § and Figure (shaded regions 
and black lines for the full and restricted dataset respectively). 

Observations are compatible with a large patch of parameter space in the a, /3/a plane 
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5 10 15 0.24 0.26 0.28 0.30 0.32 0.34 

a Q„ 



Figure 7: Likelihood Contours for disformal quintessence. Shaded regions correspond to the 1 and 
2-(T obtained by the fuU MCMC analysis while the black lines correspond to a restricted analysis 
using only information from the background expansion (SNe, BAO and CMB distance priors). The 
constraints in a and /3/a correspond respectively to high fractions of early dark energy and slow 
disformal transitions towards accelerated expansion. 



Parameter 


Maximum Likelihood 


68% C.L. 


100S7fc/i2 


2.23 


±0.05 




0.135 


±0.004 


h 


0.690 


±0.014 


a 


15.9 


>8.5 


/3/a 


17.2 


>4.9 


r 


0.086 


±0.015 




0.971 


±0.014 


Amp. 


2.91 


±0.015 



Table 2: Results from the Monte Carlo Markov Chain analysis. All the usual cosmological param- 
eters lie within one sigma from the WMAP 7 constraints presented in [|]. Amplitude is given by 
ln(10i"^) - 2r. 

and confirm the preference towards higher values of both parameters where the A limit 
can be achieved. Constraints are tighter on a due to the cumulative effect that early dark 
energy exerts in the growth of perturbations, as can be inferred by a comparison between 
the runs using the complete and the restricted sets of observations (top plot in Figure 0) . 
This difference is likely to arise when the allowed bias factor fails to compensate the effect of 
early dark energy to produce a good fit to LRG, as well as due to the effect in the third and 
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further acoustic peaks (Figures |6| and ^) . The constraints on the properties of accelerated 
expansion (i.e. on /3/a) do not improve significantly in the full analysis because the relevant 
features affect mostly distance information, encoded in SNe, BAO and CMB angular shift 
measurements. Other features such as the integrated Sachs- Wolfe effect are significant only 
for very low values of the parameters and do not introduce relevant differences in the region 
of interest (Figure ||). There is a slight mismatch between both runs in the determination 
of Om (second plot of Figure ^) . It seems unlikely that the discrepancy is caused by passing 
from the full WMAP7 spectra to the homogeneous distance priors z*,/yi,7^. Though less 
restrictive, distance priors should essentially agree with the complete data. The shift would 
be due to the additional information of the matter power spectrum from LRG. The lower 
value of pm would raise the matter power whenever the allowed linear bias is not large 
enough, since Pk oc 



4. Conclusions 

Disformal generalizations of the conformal transformation have found several applications 
in cosmology, particularly in the frameworks of gravitational alternatives to dark matter 
and varying light speed alternatives to inflation. In this paper we have shown that the 
disformal relation can be also very prolific in generating alternative explanations for the 
cosmic acceleration, suggesting new links between models worth further explorations. In 
the future, it would be interesting to study the disformal relations which may exist between 
more general field theories. Already the details of motion of point particles in a disformal 
metric are, to our knowledge, not very well understood. Whilst the nature of dark energy is 
undisclosed and it is not known to which metric it is coupled, there is in fact no fundamental 
principle establishing the precise relation between the gravitational metric and the physical 



metric in which ordinary matter lives |6^, 28, It is thus useful to investigate 

the observational signatures and put experimental bounds on it from e.g. classical tests of 
the equivalence principle. 

Here we considered the canonical scalar field action where metric is replaced by the 
disformal one, i.e. made the substitution g^jj — )■ g^y in the Lagrangian. At various limits 
we can then obtain Chaplygin gas and other models considered in the literature. Disformal 
quintessence arises in a very simple way as an application of this substitution, which results 
effectively in a non-standard self-interaction of the scalar field. This self-interaction causes 
the field to accelerate the universe given exponential functions of the field and under the 



condition a/ f3 < 1 in (^,3.2). The acceleration can begin near the present epoch if both 



a and (3 are of order one and (px ~ Mp, meaning that all the parameters appearing in the 
Lagrangian can be set to the Planck scale. For the magnitude of the potential and the 
disformal factor this follows from the special shift symmetry property of the exponential 
form. Moreover, as well known, in cosmology the exponential potential has the special 
scaling property, so practically no tuning whatsoever of initial conditions for the field is 
required. These observations motivated us to consider the phenomenology of the model in 
detail and to confront it with the latest available data. 
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We considered evolution both at the level of background and of linear perturbations. 
The dependence of the background dynamics on model parameters allows us to put con- 
straints on them. These become more stringent if we take into account the evolution of 
linear structures, though the effect of the perturbations of the field to them is more subtle. 
A Monte Carlo Markov Chain simulation was used to obtain constraints on the parameter 
space of the theory by comparing with WMAP 7 year data, SNe from the Union dataset, 
baryon acoustic peak position and luminous red galaxies power spectrum from SDSS. Six 
cosmological parameters were allowed to vary together with a and /3. 

It was then found that small values of a can be ruled out because of the effect of 
early dark energy that implies. On the other hand, small values of (3 /a result in slow 
transition to acceleration, which is also disfavored by the data. However, when both a and 
/3 /a are sufficiently large (the precise limits shown in Figure |^ , it becomes very difficult 
to distinguish the model from ACDM and thus this parameter region remains compatible 
with the present data. These constraints might be improved by considering nonlinear scales 
and gaining a deeper understanding of the effects on large scale structure. 

Further research might focus on the role of the disformal relation in dark energy sce- 
narios by the introduction of couplings to other forms of matter and gravity constructed by 
following the disformal prescription Qf^i, — )• as well the relation of disformal quintessence 
with other forms of dark energy and non- minimal derivative coupled version of scalar tensor 
theories. 
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